home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / zgbcon.z / zgbcon
Text File  |  1998-10-30  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZGGGGBBBBCCCCOOOONNNN((((3333FFFF))))                                                          ZZZZGGGGBBBBCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZGBCON - estimate the reciprocal of the condition number of a complex
  10.      general band matrix A, in either the 1-norm or the infinity-norm,
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND, WORK,
  14.                         RWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, KL, KU, LDAB, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          INTEGER        IPIV( * )
  23.  
  24.          DOUBLE         PRECISION RWORK( * )
  25.  
  26.          COMPLEX*16     AB( LDAB, * ), WORK( * )
  27.  
  28. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  29.      ZGBCON estimates the reciprocal of the condition number of a complex
  30.      general band matrix A, in either the 1-norm or the infinity-norm, using
  31.      the LU factorization computed by ZGBTRF.
  32.  
  33.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  34.      condition number is computed as
  35.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  36.  
  37.  
  38. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  39.      NORM    (input) CHARACTER*1
  40.              Specifies whether the 1-norm condition number or the infinity-
  41.              norm condition number is required:
  42.              = '1' or 'O':  1-norm;
  43.              = 'I':         Infinity-norm.
  44.  
  45.      N       (input) INTEGER
  46.              The order of the matrix A.  N >= 0.
  47.  
  48.      KL      (input) INTEGER
  49.              The number of subdiagonals within the band of A.  KL >= 0.
  50.  
  51.      KU      (input) INTEGER
  52.              The number of superdiagonals within the band of A.  KU >= 0.
  53.  
  54.      AB      (input) COMPLEX*16 array, dimension (LDAB,N)
  55.              Details of the LU factorization of the band matrix A, as computed
  56.              by ZGBTRF.  U is stored as an upper triangular band matrix with
  57.              KL+KU superdiagonals in rows 1 to KL+KU+1, and the multipliers
  58.              used during the factorization are stored in rows KL+KU+2 to
  59.              2*KL+KU+1.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZGGGGBBBBCCCCOOOONNNN((((3333FFFF))))                                                          ZZZZGGGGBBBBCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      LDAB    (input) INTEGER
  75.              The leading dimension of the array AB.  LDAB >= 2*KL+KU+1.
  76.  
  77.      IPIV    (input) INTEGER array, dimension (N)
  78.              The pivot indices; for 1 <= i <= N, row i of the matrix was
  79.              interchanged with row IPIV(i).
  80.  
  81.      ANORM   (input) DOUBLE PRECISION
  82.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  83.              NORM = 'I', the infinity-norm of the original matrix A.
  84.  
  85.      RCOND   (output) DOUBLE PRECISION
  86.              The reciprocal of the condition number of the matrix A, computed
  87.              as RCOND = 1/(norm(A) * norm(inv(A))).
  88.  
  89.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  90.  
  91.      RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
  92.  
  93.      INFO    (output) INTEGER
  94.              = 0:  successful exit
  95.              < 0: if INFO = -i, the i-th argument had an illegal value
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.